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Abstract. We have modelled radio and far-infrared rotational 
lines of 12 CO and 13 CO obtained from the circumstellar en- 
velope of the infrared-bright carbon star IRAS 15194-51 15. 
Eleven rotational lines between J=l-0 and 21-20 and nine 
rotational lines between J=l-0 and 22-21 in the ground vi- 
brational states of 12 CO and 13 CO, respectively, provide the 
observational constraints. A model of the circumstellar enve- 
lope with a constant mass-loss rate (1 x 10~ 5 M Q yr _1 ) and 
12 CO/ 13 CO-ratio (5.5) is consistent with our observed FIR and 
radio data, suggesting that the wind characteristics have not 
changed significantly over the past few thousand years. Thus, 
IRAS 15194-5115 appears to be a highly evolved AGB-star, 
but the carbon star properties combined with the inferred low 
12 C/ 13 C-ratio make the evolutionary status of this star uncer- 
tain. It may have been a J-star for which the 12 C/ 13 C-ratio has 
remained low, or it may be a star of 5 to 8 solar masses, which 
has recently become a carbon star due to quenching of hot bot- 
tom burning. The dust properties or the relative amount of dust 
in the envelope appear to be different from those in the enve- 
lope of the well known carbon star IRC+10 216. 

Key words: stars: AGB and post-AGB - circumstellar matter 
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1. Introduction 

Red giant stars lose large amounts of matter through a mas- 
sive, but slow, stellar wind during their late evolution on the 
Asymptotic Giant Branch (AGB). The mass loss results in cir- 
cumstellar envelopes (CSEs), several thousand stellar radii in 
size, consisting of gas, mainly in molecular form, and dust. Im- 
portant information on late stellar evolution can be obtained 
from the circumstellar emission (Olofsson 1996). 

We have studied the third brightest (at 12 /im) carbon star 
(C-star) IRAS 15194-5115 (or II Lup), which is the brightest 
one known in the southern hemisphere, in several rotational 
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transitions in the ground vibrational states of 12 CO and 13 CO. 
The lowest energy lines, J=l-0 etc., lie at radio wavelengths 
and probe the outer regions of the CSE. Higher energy lines, 
for example, the J=14-13 to J=22-21 lines presented in this 
paper, which lie in the far-infrared (FIR) wavelength range 
100-200 /im originate from regions much closer to the star. 
By combining observations in these two wavelength ranges we 
are able to probe different regions of the CSE, and to compare 
wind characteristics, such as mass loss rate, expansion velocity, 
molecular abundances (including isotope ratios), over a time 
span of several thousand years. 

IRAS 15194-5115 was discovered by Epchtein et al. (1987; 
the ESO-Valinhos 2.2 /im sky survey) and by IRAS in the far- 
infrared, and it was identified as a C-star using optical and in- 
frared spectra (Meadows et al. 1987). To understand the late 
evolutionary phases of C-stars, in which they undergo intense 
mass loss, this star provides an object which is bright enough to 
make possible an exhaustive spectroscopic investigation over 
a broad range of wavelengths. IRAS 15194-5115 shows sim- 
ilarities in the infrared to the well-studied prototype C-star 
IRC+10216 (Meadows et al. 1987; Lopez et al 1993). On 
the basis of radio observations Nyman et al. (1993) concluded 
that most molecular species, except a few, have circumstellar 
abundances comparable to those in IRC+10216. Also, they 
found that the two objects have comparable mass loss rates. 
The most striking difference between these two C-stars is the 
12 C/ 13 C-ratio, which appears to be at least five times lower in 
IRAS 15 194-5 115 (Nyman et al. 1993). 

Here, we present observations of FIR rotational lines of 
12 CO and 13 CO from IRAS 15194-5115, as well as a model 
interpreting these data and radio rotational lines. 

2. Observations 

We have obtained FIR observations with the Infrared Space 
Observatory (Kessler et al. 1996), see Fig. |IJ The Long 
Wavelength Spectrometer (LWS, Clegg et al. 1996), with a 
beam size assumed to be 70", was used in the grating mode 
(LWS01), covering the range 90-197 /im, and providing a mean 
spectral resolution element of ^0.7 /im (this corresponds to 
R = A/AA ~ 100 — 250; hence the circumstellar lines are far 
from being resolved). The spectra were sampled at 0.15 /im. 
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Fig. 1. The FIR spectrum of IRAS 15194-51 15 as observed by the ISO satellite. Fluxes from four detectors are shown without 
scaling. The overlap regions show the uncertainties, both the relative and the absolute ones, in these ISO-LWS observations. All 
rotational transitions within the ground vibrational states of 12 CO, 13 CO, H 12 CN, and H 13 CN in this region are indicated. Also 
the interstellar [Cll] line is marked. These emission lines are superimposed on the continuum emitted by the dust around the star 



The 3.5 h of observations were performed on 19 September 
1997 during ISO revolution 674. The reductions were made 
using the most recent pipeline basic reduction package OLP 
(v.7) and the ISO Spectral Analysis Package (ISAP v. 1.5). The 
pipeline processing of the data, such as wa veleng th and flux 
calibration, is described in Swinyard et al. ( 1996 ); the com- 
bined absolute and systematic uncertainties in the fluxes are of 
the order of ±30%. The accuracy of the wavelength varies and 
can be as bad as ±0. 1 /im. The observed fluxes are measured 
after scaling the four detectors and averaging the two measure- 
ments in the overlap regions between different LWS detectors. 
The scaling needed is at most a few percent. In Fig. 1 the fluxes 
from the four detectors are shown without any scaling. Addi- 
tional uncertainties in the line fluxes, sometimes large, are in- 
troduced in the continuum subtraction. A dust continuum level 
is estimated and subtracted over a limited wavelength interval, 
usually ~10/im. The errors in the fluxes due to the uncertain 
continuum subtraction are estimated to be around ±20%. A 
few CO-lines are blended by an unknown amount of H 13 CN, 
leading to somewhat too high measured CO-fluxes in these 
lines. 

Our final LWS spectrum is dominated by thermal radiation 
due to the large quantity of dust around the o bject. The dust 
modelling of the object by Lopez et al. ( 1993 ) shows that for 
wavelengths longer than 10 /im the spectrum is totally dom- 



inated by dust emission. Superimposed on this emission are 
molecular emission lines, of which we have identified those due 
to 12 CO, 13 CO, and H 12 CN (all relevant rotational transitions 
within the ground vibrational states of 12 CO, 13 CO, H 12 CN, 
and H 13 CN are indicated in Fig. [[]). Here we present eight 
lines between J=14-13 to 21-20 of 12 CO and seven lines be- 
tween J=15-14 and 22-21 of 13 CO, Table 1 and Fig. |} A sim- 
ilar LWS spectrum of IRC+10216 has been presented by Cer- 
nicharo et al. (1996). The observed fluxes from this source are 
an order of magnitude higher than those towards IRAS 15194- 
5115, making the analysis more difficult in our case. 

As a supplement to the ISO data, we have used 12 CO and 
13 CO radio-line observations from SEST [Swedish-ESO Sub- 
millimetre Telescope; the J=l-0 and 2-1 spectra published in 
Nyman et al. (1993), and a 12 CO(3-2) spectrum (Nyman, pri- 
vate communication)]. The radio data are presented in Fig. ^. 
The errors in the radio intensities are about ±10%, and the res- 
olution is better than 1 km/s. 

3. The model 

In order to model the circumstellar emission lines we have 
used a numerical simulation code, based on the Monte Carlo 
method, to solve the non-LTE radiative transfer of 12 CO and 
13 CO taking into account 30 rotational levels for each of the 
ground and first vibrational states (see Schoier, PhD thesis in 
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Table 1. Observed and modelled CO FIR rotational lines from 
IRAS 15194-5115. A colon (:) marks lines with uncertain flux 
estimates due to a complex continuum around the line around 
143.5/^m or due to a blending of two or three lines of 12 CO, 
13 CO or HCN 
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prep., for details). Using 40 energy levels does not change the 
temperature structure and the line intensities noticeably. For 
example, in the innermost regions, cooling by CO excited to 
higher levels than J=20-25 are unimportant. Crosas & Menten 
(1997) also found that around 30 levels are sufficient to treat 
the CO cooling in CSEs properly. The decline in CO cooling in 
the inner regions of the CSE, as demonstrated in Fig. 4, is due 
to the fact that the CO lines become optically thick. 

The transition p robabi lities and energy levels are taken 
from Chandra et al. ( |l996| ), and the rotational (C O-H 2 ) colli- 
sional rate coefficients are from Flower & Launay ( 1985 ) (they 
are extrapolated for J>11 and for temperatures higher than 
250 K). The derived intensities of the FIR lines are not very 
sensitive to the accuracy of this extrapolation since they are 
formed in LTE (see below). Collisional transitions between vi- 
brational levels are not important due to the low density and the 
short radiative lifetimes. 

The kinetic gas temperature is calculated in a self- 
consistent way, using the derived level populations, by solving 
the energy balance equation (e.g. Goldreich & Scoville 1976), 



dT 
dr 



(2 - 2 7 ) 



T 



7 



1 



(H-C), 



(1) 



where 7 is the adiabatic index, ?t,h 2 the number density of 
molecular hydrogen which is the most abundant species in the 
CSE, Vqc the terminal gas expansion velocity, k the Boltzmann 
constant, H the total heating rate per unit volume, and C the 
total cooling rate per unit volume. 

The first term on the right-hand-side in the energy balance 
equation is the cooling due to the adiabatic expansion of the 
gas. The adiabatic index, i.e. the ratio of the specific heats 
7 = c p /c v , can be assumed to be 5/3 since H 2 can be treated 



as a 'mono-atomic gas'. This is true as long as the corrections 
on 7 due to rotational excitation of H2 are small. These correc- 
tions are important for temperatures greater than approximately 
300 K, i.e. for the innermost regions where 7=7/5 (Goldreich 
& Scoville 1976). According to our tests, the effect of a varied 
adiabatic index (in comparison to a constant one) is only a few 
percent in the intensities, which is in agreement with the result 
of Groenewegen (1994). 

Additional cooling is provided by emission from the 
molecules CO and H2. The molecular cooling due to CO ( 12 CO 
and 13 CO) is calculated from the derive d level populations us- 
ing the expression of Crosas & Menten (1997). For H2 we use 
the approach by Groenewegen (1994). HCN could be an im- 
portant coolant in the inner parts of the envelope (Cernicharo 
et al. 1996), but Groenewegen (|1994|) demonstrates that HCN 



cooling is only of minor importance. In the present version of 
the code HCN cooling is therefore not included. 

The heating is assumed to be dominated by dust-gas 
collisions (Groenewegen 1994) and the photoelectric effect 
(Crosas & Menten 1997). In the former case a number of 
free parameters describing the dust are introduced. These 
are highly uncertain but affect the derived line intensities. 
Here we assume that the Q-parameter, i.e., the efficiency 
of momentum transfer, to be equal to 0.03, and define a 
new parameter that contains the other dust parameters, h = 
(*/0.01)(2.0gcm"" 3 / / 9 d )(0.05^m/a d ), where ^ is the gas- 
to-dust ratio, pd is the dust density, and the mean size of the 
dust particles. The normalised values are the ones we used to 
fit the CO line emission of IRC+10216 using our model (see 
Schoier et al. in prep.). Heating (or cooling) due to a temper- 
ature difference (heat exchange) between the dust and the gas 
could affect the derived temperature structure of the gas and 
hence the line intensities, especially in the inner parts where 
the coupling between dust and gas is high. However, using the 
dust temperature structure for IRAS 15194-5115 from Lopez 



et al. (1993) and the analytical expression from Groenewegen 
(1994), we conclude that this effect is only of minor impor- 
tance. 

IRAS 15194-5115 is a LPV with a period of 580 days, and 
with a variation of 0.8 mag in the Johnson M filter (Le Bertre, 
1992), but the amplitude decreases rapidly towards longer 
wavelengths. The period-luminosity relation for C-stars (Groe- 
newegen & Whitelock 1996) gives a luminosity of 8770 L . 
Using this luminosity and the apparent bolometric flux we ar- 
rive at a distance of 600 pc, a value in agreement with the ones 
discussed recently in the literature. The stellar temperature is 
set to T c ff = 2200 K, a value typical for C-stars with a high 
mass loss rate (e.g., Le Bertre 1997). However, the mass loss 
rate derived is quite insensitive to the choice of T e g, since ra- 
diative excitation is of very little importance in this case. The 
reason is that, as our calculations show, in this high-density 
CSE the populations in the CO rotational levels are deter- 
mined by collisions, and hence the kinetic temperature, out to 
a fewx 10 16 cm (Fig. Q). Consequently, the light variability of 
the object is also unimportant for the excitation. In the present 
version of our code only the star provides radiation since no 
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Fig. 2. Millimetre-wave CO spectra (histogram) observed towards IRAS 15194-5115 using SEST, overlaid with the best model 
results (full line) using a mass loss rate of 1 x 10~ 5 M yr _1 and a 12 CO/ 13 CO ratio of 5.5. The beam size of each observation 
is given 



thermal emission from dust in the CSE is included. This latter 
radiation could be significant for pumping vibrational states, 
leading to non-thermal level populations of the rotational states 
which could effect the FIR line intensities, but this effect is 
judged to be small in this high-density CSE. 

We assume a spherically symmetric CSE formed by a con- 
stant mass loss rate and expanding at a constant velocity. From 
the widths of the radio lines we are able to estimate the ter- 
minal gas expansion velocity, ~21.5kms _1 . The turbulent ve- 
locity width affects the effectiveness of the radiative excitation. 
However, since this is only of minor importance, we will use a 
constant value of 0.5 km s -1 throughout the envelope. The in- 
ner radius we set to 4 R* ; the stellar radius being estimated to 
be 4.5xl0 13 cm [we note here that Lopez et al. (1993) found 
the dust condensation radius to be about 6R*, where the dust 
temperature is ~1000 K]. The inner parts of the CSE are com- 
plicated regions of which we have a very limited understand- 
ing, see e.g. Ryde et al. ( 1999| ), and where many of the basic 
assumptions of our model are certainly not valid. The prop- 
erties of this region could, in principle, be of importance for 
the excitation of the higher rotational levels. According to our 
tests, however, an initial velocity gradient for example, affects 
the intensities only marginally. We have also performed tests 



showing that a change of the inner radius by a factor of two af- 
fects the derived FIR line intensities by at most a few percent. 

The photodissociation radius of the CO envelope is defined 
as the location where the fractional abundance of CO is half the 
photospheric value. This was estimated to be 3.2 x 10 17 cm for 
12 CO (for the esti mated mass loss rate) using the calculations 
by Mamon et al. ( 1988 ). The two isotopic variants have dif- 
ferent photodissociation and c hemica l fractionation behaviour, 
but according to Mamon et al. ( 1988 ) these effects tend to can- 
cel out, and the authors predict that the 12 CO/ 13 CO-ratio is 
constant with radius, except that the 13 CO envelope extends 
slightly further out. Here we assume the same radius for both 
isotopic variants. The J=l-0 lines are affected the most when 
varying the outer radius, but according to our calculations only 
by ~20% in intensity, even if the outer radius is changed by a 
factor of two. We adopt an initial abundance of 12 CO with re- 
spect to H2, /co, of 10~ 3 . The estimated mass loss rate scales 
roughly as Jqq 8 - The 13 CO abundance is obtained from the 
best-fit model. 



4. Results and discussion 

When running our model we vary the free parameters, the mass 
loss rate, the 12 CO/ 13 CO-ratio, and the dust parameter h, until 
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Fig. 3. Selected continuum-subtracted FIR CO spectra (histogram) observed towards IRAS 15194-5115 using the ISO satellite, 
overlaid with the best model results (full line) using a mass loss rate of 1 x 10 -5 M Q yr" 1 and a 12 CO/ 13 CO ratio of 5.5. For all 
transitions the beam size and spectral resolution used are 70" and 0.7/im, respectively 



we get a good fit to the observed 12 CO and 13 CO radio line 
intensities, profiles, and radial brightness distributions (Nyman 
et al. 1993). This gives the kinetic temperature structure of the 
CSE, albeit constrained only by emission from the outer re- 
gions. Finally, we compare the model fluxes for the higher ro- 
tational transitions with our LWS observations. In Figs [| & |] 
we have plotted the line profiles (convolved with the instrumen- 
tal resolutions) of the best-fit model to the 12 CO lines (there are 
narrow interstellar CO lines present in the blue wing of the ra- 
dio lines, and we exclude these lines when integrating the radio 



fluxes). We also plot the best model fit to the 13 CO lines for the 
best estimate of the 12 CO/ 13 CO-ratio. 



From a simple error analysis, we arrive at a mass loss rate 
(based on the radio observations) of (l±0.2)x 1O -5 M0 yr _1 
(excluding uncertainties in distance, 12 CO abundance, and the 
assumed circumstellar model), in agreement with other es- 
timates (Nyman et al. 1993; Le Bertre 1997). We obtain a 
12 CO/ 13 CO-ratio of 5.5±1. This result agrees excellently with 
that found by Nyman et al. ( 1993 ) using intensity ratios of iso- 
topic variants of CS, CN, HCN, and HC 3 N. 
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Fig. 4. The full line in the temperature panel gives the gas kinetic temperature, the dashed [dotted] line is the excitation temper- 
ature of the 12 CO(1-0) [ 12 CO(18-17)] transition. The dash-dot [dash-dot-dot-dot] line gives the tangential optical depth, r tan , 
of the 12 CO(1-0) [ 12 CO(18-17)] transition. In the cooling panel the full line represents adiabatic cooling; the dotted line the H2 
cooling; and the dashed line and the dash-dot line the 12 CO and 13 CO cooling, respectively (see text for details) 



The kinetic and excitation (for selected transitions) temper- 
ature structures are plotted in Fig. Q The 12 CO(1-0) transi- 
tion is close to thermalised out to ^3xl0 16 cm. The 13 CO(l- 
0) transition is strongly super-thermal in the inner parts, and it 
is never thermalised. The FIR lines are thermalised only out 
to about lxl0 15 cm. The radio lines are formed mainly be- 
yond ~5xl0 16 cm (^1000R*), corresponding to a time scale 
of JblO 3 yr, while the FIR lines are formed much closer to the 
star, within 10 15 cm (~10R*), which corresponds to ^20 yr. 

In Fig. fi| we also plot the radial variation of the various 
cooling terms included in the energy balance equation, Eq. ([!]). 
Here, for clarity, the cooling rate has been multiplied by a fac- 
tor oc r 4 (the heating due to dust-gas collisions goes roughly 
as r~ 4 ). CO cooling (mainly by 12 CO) is the dominant coolant 
in the region between ^4xl0 14 cm and ^lxl0 17 cm, where 
most of the emission from the observed transitions emanates. 
In this context a high abundance of 13 CO could affect the 
modelled line intensities significantly. For IRAS 15194-5115 
a 12 CO/ 13 CO-ratio of only 5.5 implies that cooling by 13 CO 
could be as high as 20% of that of 12 CO. We have inluded 
13 CO cooling in the same self-consistent way as for 12 CO and 
it turns out that its relative cooling is somewhat higher than 
20%, Fig. |, due to the fact that the lines of 13 CO are opti- 
cally thinner than those of 12 CO. Closer to the star H 2 cooling 
dominates but very little emission from this part of the enve- 
lope contributes to the line emission observed. In the cool, ten- 
uous, outer parts of the envelope, where most of the observed 
CO(l-O) emission arises, cooling due to the adiabatic expan- 
sion dominates. 

We derive a dust parameter h^2 (with an estimated error 
of about 30%), i.e. twice as large as the one we used to fit the 



CO radio lines of IRC+10216 in our modelling of this object. 
This indicates that the dust properties or the relative amount of 
dust could be different in the two objects. Our result is consis- 
tent with the finding by Groenewegen et al. (1998 ) in their dust 
models that IRAS 15194-5115 has a dust-to-gas ratio which is 
about twice the value for IRC+10216 (using the same values 
for other dust parameters). However, our derived dust parame- 
ter is sensitive to the adopted luminosity, and hence the result 
is uncertain. 

The model fluxes compare well with the observed FIR 
fluxes within the estimated errors, see Table 1, although we 
note that the fluxes from lines at longer wavelengths are sys- 
tematically higher compared to what the model predicts, pos- 
sibly indicating a systematic error. We conclude that the FIR 
data, probing the inner regions of the CSE, are consistent with 
the model derived from radio observations. Thus, we suggest 
that there have been no significant changes in the stellar wind 
characteristics over the past ^5000 years. Due to the estimated 
large uncertainties in the LWS data the mass loss rate in the in- 
ner parts may have changed by about a factor of three and still 
be consistent with the model. Likewise, there appears to have 
been no significant change in the 12 CO/ 13 CO-ratio over this 
period. 

We believe that the circumstellar 12 CO/ 13 CO-ratio is 
very close to the 12 C/ 13 C-ratio in the stellar atmosphere. A 
12 C/ 13 C-ratio of only 5.5 is just slightly larger than the equi- 
librium value of 3-3.5 found from hydrogen burning through 
the CNO-cycle. Such a low 12 C/ 13 C-ratio puts IRAS 15194- 
5115 close to the group of C-stars known as the J-type stars, 
which have 12 C/ 13 C-ratios of 3-4 (Lambert et al. 1986). Opti- 
cally bright C-stars often have ratios above 30 (Lambert et al., 
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1986), and this seems to apply also for high mass loss rate ob- 
jects (Greaves & Holland 1997). The best example of the lat- 



ter is IRC+10216 for which Kahane et al. ( |1992[ ) derived a 
value of 12 C/ 13 C=44±3. From its 12 C/ 13 C-ratio IRAS 15194- 
5115 resembles more closely the extreme OH/IR-stars which 
are massive AGB stars in the very last phases of evolution 
(Delfosse et al. 1997). The low 12 C/ 13 C-ratio of the OH/IR 
stars may be caused by hot bottom burning. This occurs when 
the temperature is larger than ^80xl0 6 K and this is only 
achieved in massive (5-8 M Q ) AGB stars (Boothroyd et al., 
1993). The luminosity derived for IRAS 15194-5115 from the 
period-luminosity relation (8770 L ) is, however, lower than 
typical OH/IR luminosities (-3 x 10 4 L ). Thus, IRAS 15194- 
5115 may have been a J-star for which the 12 C/ 13 C-ratio has 
remained low despite it probably being highly evolved. Alter- 
natively, it may be a fairly massive star, which has recently be- 
come a carbon star due to quenching of hot bottom burning 
(Frost et al. 1998). 

We have shown that modelling radio and FIR lines simul- 
taneously is potentially an interesting diagnostic method for 
probing the stellar wind history of giant stars. Of special in- 
terest are, e.g., the group of OH/IR stars which are suggested 
to be undergoing an onset of a superwind phase with a drasti- 
cally enhanced mass-loss rate during the recent 500-1000 years 
(Delfosse et al. 1997). This hypothesis could be observationally 
confirmed either with interferometric CO observations, as sug- 
gested by Delfosse et al. (1 997 ) or by FIR observations of lines 



from rotationally excited CO molecules as is demonstrated in 
this paper. 
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